Molecular Plant

Fruit ripening evolved to be attractive to frugivores that derive energy and nutrition from the fruits in exchange for assisting seed dispersal, which is accompanied by the dramatically change of fruit characteristics, including color, aroma, and texture. The plant hormone ethylene plays a key role in climacteric fruit ripening, while the role of other phytohormones as well as their cross talk with ethylene in modulating fruit ripening remains elusive. Here, we report growth-promoting phytohormone brassinosteroids promote fruit ripening in tomato through regulation of ethylene biosynthesis. Exogenous BR treatment and the increase of endogenous BR content in SlCYP90B3-OE promoted ethylene production and fruit ripening. SlBZR1, a central component and positive regulator of BR signaling pathway, promotes ethylene production and carotenoid accumulation through direct transcriptional regulation of SlACO1, SlACO3 and SlPSY1. Furthermore, SlBIN2, a negative regulator of BR signaling upstream of SlBZR1, decreases ethylene production and carotenoid accumulation. Together, our results demonstrate that BR signaling integrates ethylene and carotenoid biosynthetic pathway to regulate fruit ripening.

Sessile plants exhibit diverse movement behaviors that have long intrigued the scientists. The legume plants display a rhythmic leaflet movement pattern characterized by horizontal opening during the day and vertical closure at night. However, the underlying mechanisms remain largely enigmatic. Here, we isolated a mutant designated as dlm1 (downward leaflet movement1) from Medicago truncatula that displays leaflets downward opening during daytime while upward closure at night. Cellular analyses reveal that this aberrant phenotype correlates with abnormal volume changes of motor cells within the pulvinus. DLM1 exhibits high expression level in the motor organ pulvinus and encodes an unreported nuclear-localized PP2C phosphatase. Notably, phylogenetic analysis demonstrates that species exhibiting rhythmic leaflet movements consistently retain DLM1 homologs, suggesting the possibility of its functionally conserved role in regulating leaflets movement pattern. Structural characterization reveals that DLM1 possesses both phosphatase and kinase domains. Functional complementation assays demonstrate that the phosphatase domain is necessary and sufficient for maintaining the leaflet movement pattern. Collectively, our work uncovers a novel PP2C protein that governs the leaflet movement, providing mechanistic insights into this intriguing plant behavior.

The Asteraceae (Compositae) is the largest flowering plant family, ubiquitous in most terrestrial communities, and morphologically hyper-diverse. An ancient whole genome triplication (paleo-hexaploidization) occurred at approximately the same time as the evolutionary innovation and adaptive radiation of the family during the middle Eocene. Despite its importance, the genomic contents arising from this triplication have yet to be tracked in context of the Asteraceae genome evolution. We applied a synteny oriented phylogenomic analysis of 21 Asterales genomes and to study the paleo-hexaploidization and its consequences to gene, trait, and genome evolution. We identified 15 ancestral linkage groups (ALGs) that date back to the common diploid ancestor of all Asteraceae. Each of these groups was triplicated, resulting in 45 genomic blocks (3x15), which serve as the foundation for cross-family analyses. We demonstrate the complex evolutionary dynamics of the 45 genomic blocks across the Asteraceae phylogeny. We found that modern genomes are genetic mosaics of three progenitor genomes by extensive genomic exchange, chromosomal shuffling and gene fractionation. 157 genes retained three paleo-hexaploid derived syntenic paralogs across most Asteraceae species. Transcription factors (TFs) and auxin-related genes are significantly overrepresented in the conserved triplets, and expression of the paleo-hexaploidy paralogs is spatiotemporally differentiated. These genes are involved in the development of floral capitulum, a remarkable morphological innovation of the family. The discovery of conserved triplicated genes can direct further study to understand the evolutionary innovation, and the synteny-phylogenomic framework and ALGs provide a comparative framework to characterize newly sequenced Asteraceae genomes.

Cross-regulation between hormone signaling pathways is indispensable for plant growth and development. However, the molecular mechanisms by which multiple hormones interact and co-ordinate activity need to be understood. Here, we generated a cross-regulation network explaining how hormone signals are integrated from multiple pathways in etiolated Arabidopsis (Arabidopsis thaliana) seedlings. To do so we comprehensively characterized transcription factor activity during plant hormone responses and reconstructed dynamic transcriptional regulatory models for six hormones; abscisic acid, brassinosteroid, ethylene, jasmonic acid, salicylic acid and strigolactone/karrikin. These models incorporated target data for hundreds of transcription factors and thousands of protein-protein interactions. Each hormone recruited different combinations of transcription factors, a subset of which were shared between hormones. Hub target genes existed within hormone transcriptional networks, exhibiting transcription factor activity themselves. In addition, a group of MITOGEN-ACTIVATED PROTEIN KINASES (MPKs) were identified as potential key points of cross-regulation between multiple hormones. Accordingly, the loss of function of one of these (MPK6) disrupted the global proteome, phosphoproteome and transcriptome during hormone responses. Lastly, we determined that all hormones drive substantial alternative splicing that has distinct effects on the transcriptome compared with differential gene expression, acting in early hormone responses. These results provide a comprehensive understanding of the common features of plant transcriptional regulatory pathways and how cross-regulation between hormones acts upon gene expression.

Momilactone A, an important plant labdane-related diterpenoid, functions as a phytoalexin against pathogens and an allelochemical against neighboring plants. The genes involved in biosynthesis of momilactone A are found in clusters, i.e., MABGCs (Momilactone A biosynthetic gene clusters), in the rice and barnyardgrass genomes. How MABGCs originate and evolve is still not clear. Here, we integrated results from comprehensive phylogeny and comparative genomic analyses of the core genes of MABGC-like clusters and MABGCs in 40 monocot plant genomes, providing convincing evidence for the birth and evolution of MABGCs in grass species. The MABGCs found in the PACMAD clade of the core grass lineage (including Panicoideae and Chloridoideae) originated from a MABGC-like cluster in Triticeae (BOP clade) via horizontal gene transfer (HGT) and followed by recruitment of MAS and CYP76L1 genes. The MABGCs in Oryzoideae originated from PACMAD through another HGT event and lost CYP76L1 afterwards. The Oryza MABGC and another Oryza diterpenoid cluster c2BGC are two distinct clusters, with the latter being originated from gene duplication and relocation within Oryzoideae. Further comparison of the expression patterns of the MABGC genes between rice and barnyardgrass in response to pathogen infection and allelopathy provides novel insights into the functional innovation of MABGCs in plants. Our results demonstrate HGT-mediated origination of MABGCs in grass and shed lights into the evolutionary innovation and optimization of plant biosynthetic pathways.

Plants rely on the innate immune system to sense and respond to a wide range of lifestyle pathogens and to facilitate their survival in natural ecosystems. Pathogen-associated molecular patterns (PAMP)-triggered immunity (PTI) and effector-triggered immunity (ETI) are designated as a two-branched system of innate immunity. Although PTI/ETI share a series of downstream molecular events, systematic analysis of convergent and divergent signaling in PTI/ETI is currently lacking. The phytohormones salicylic acid (SA) and jasmonic acid (JA) are considered to constitute the hormonal backbone of plant immunity, are functionally antagonistic, and play essential roles in defending against biotrophic and necrotrophic pathogens, respectively. However, the distinct performance of two phytohormones in PTI/ETI remains unclear. Here, we systemically investigate and validate the reprogramming of molecular networks during PTI and ETI. Using publicly available Arabidopsis RNA sequence data from 560 samples, we construct a co-expression network under Mock conditions and then explore the differential expression/co-expression changes during PTI, ETI, and Pto DC3000 infection. During PTI and ETI, one-third of genes in the Arabidopsis genome exhibit the same directional differential expression in a manner independent of JA/ethylene/PAD4/SA signaling but show differential co-expression patterns. However, the defense phytohormone network is required for defense against Pto DC3000 infection. We also exhibit the use of this network in prioritizing genes that functioned closely with the proteins directly targeted by elicitors. Overall, this study will deepen our understanding of plant transcriptome in plant immunity and provide new insights into the mode of action of elicitors.

GOLDEN2-LIKE (GLK) transcription factors drive the expression of photosynthesis-associated nuclear genes (PhANGs), indispensable for chloroplast biogenesis. We previously demonstrated that the salicylic acid (SA)-induced SIGMA FACTOR-BINDING PROTEIN 1 (SIB1), a transcription coregulator and positive regulator of cell death, interacts with GLK1 and GLK2 to reinforce their activities. The SIB1-GLK interaction raises the level of light-harvesting antenna proteins in photosystem II, aggravating photoinhibition and singlet oxygen (1O2) burst. 1O2 then contributes to SA-induced cell death via EXECUTER 1 (EX1, 1O2 sensor protein)-mediated retrograde signaling upon reaching a critical level. We now reveal that LESION-SIMULATING DISEASE 1 (LSD1), a transcription coregulator and negative regulator of SA-primed cell death, interacts with GLK1/2 to repress their activities. Consistently, the overexpression of LSD1 represses GLK target genes including PhANGs, whereas the loss of LSD1 increases their expression. Remarkably, LSD1 overexpression inhibits chloroplast biogenesis, resembling the characteristic glk1glk2 double mutant phenotype. The subsequent chromatin immunoprecipitation analysis coupled with quantitative PCR further reveals that LSD1 inhibits the DNA-binding activity of GLK1 towards its target promoters. The SA-induced nuclear-targeted SIB1 appears to counteractively interact with GLK1/2, leading to the activation of EX1-mediated 1O2 signaling. Taken together, we provide a working model that SIB1 and LSD1, mutually exclusive SA-signaling components, antagonistically regulate GLK1/2 to fine-tune the expression of PhANGs, thereby modulating 1O2 homeostasis and related stress responses.

The phytohormone auxin, specifically indole-3-acetic acid (IAA) plays a prominent role in plant development. Cellular auxin concentration is coordinately regulated by auxin synthesis, transport, and inactivation to maintain auxin homeostasis; however, the physiological contribution of auxin inactivation to auxin homeostasis has remained elusive. The GH3 genes encode auxin amino acid conjugating enzymes that perform a central role in auxin inactivation. The chemical inhibition of GH3s in planta is challenging because the inhibition of GH3 enzymes leads to IAA overaccumulation that rapidly induces GH3 expression. Here, we developed a potent GH3 inhibitor, designated as kakeimide (KKI), that selectively targets auxin-conjugating GH3s. Chemical knockdown of the auxin inactivation pathway demonstrates that auxin turnover is very rapid (about 10 min), indicating auxin biosynthesis and inactivation dynamically regulate auxin homeostasis.

Riptortus pedestris (Fabricius) a major soybean pest migrates into soybean fields during pod filling stage resulting in a leaf and stem staygreen while pods without beans syndrome. Given the agricultural importance of this species and the lack of characterized HAMP from piercing-sucking insects we performed a large scale of screening by expression of 87 R. pedestris salivary proteins with signal peptides in Nicotiana benthamiana obtaining a candidate HAMP RPH1. RPH1 activated a series of PTI responses including ROS burst upregulation of defense marker genes such as PR1 WRKY7 WRKY8 Acre31 and CYP71D20 MAPK activation and biosynthesis of phytohormones in plants. RPH1 significantly enhances soybean resistance against R. pedestris feeding. PRR coreceptors BAK1 and SOBIR1 were required for RPH1-induced PTI responses. Remarkably RPH1 homologs were widely distributed in herbivorous insects and majority of homologs from selected species induced cell death or ROS. Thus our results demonstrated that RPH1 is a conserved HAMP within chewing and piercing-sucking insects. We also discovered that R. pedestris evolved four paralogs to overcome the plant immunity triggered by RPH1. This study filled a major gap of HAMP identification from piercing-sucking insect and also deciphered a novel evasion strategy of plant immunity exploited by herbivorous insects. One sentence summaryRiptortus pedestris RPH1, a conserved HAMP in herbivores, activates a variety of PTI responses in plants. To couterdefense, R. pedestris evolved four paralogs to suppress RPH1-induced PTI responses. The author(s) responsible for distribution of materials integral to the findings presented in this article in accordance with the policy described in the Instructions for Authors (https://academic.oup.com/plcell/pages/General-Instructions) is: Ai Xia (xia@njau.edu.edu).

Receptor-like kinases (RLKs) are essential in plants and phosphorylation is a critical step for their function. Interestingly, RLKs have many non-catalytic kinases/ pseudokinases and the biochemical basis for these pseudokinases remains unclear. FERONIA (FER) is an RLK with kinase activity, but the necessity of its kinase activity for genetic functions has been debated. Here, we uncovered that the kinase-deficient variant FERK565R can activate kinase activity in FER and its homologous through homo/heterodimerization-dependent allosteric activation. We further showed that reactive oxygen species (ROS) significantly promote the dimerization of FER family members. Next, we revealed that mutating the FER P740 within the G-H loop reduces FER dimerization and disrupts its allosteric activation, thus attenuating FERs transphosphorylation for its substrate. This disruption in allosteric activation abolishes the genetic function of FERK565R, impacting ROS production and ABA-mediated stomatal movements. Additionally, we found that MEDOS1 (MDS1), a member of the FER family, is incapable of catalyzing phosphotransfer, but can boost the kinase activity of FER and HERK1 through allosteric activation. These findings settle the debate on FERs inactivated forms, and reveal a new mechanism for allosteric activation of RLKs via redox signaling, enhancing our understanding of pseudokinases in plants. One-sentence summaryFER activates kinase activity of homologous family members through allosteric activation.

Dear Editor, Since the first plant receptor-like kinase (RLK) gene ZmPK1 was cloned from Zea mays in 1990 (Walker & Zhang, 1990), this large gene family has been extensively studied and shown to play crucial roles in growth, development, and immunity (Tang et al., 2017). RLKs are widespread in the plant kingdom, while the biological functions of most RLKs remain largely elusive (Dievart et al., 2020). Given RLKs share a conserved monophyletic RLK/Pelle kinase domain, RLKs in several model plants are classified into distinct families by extracellular domains (ECDs) (Shiu & Bleecker, 2001). However, independent domain shuffling in specific lineages drives the origin of novel families, which raises a question: how about the landscape of RLKs in the whole plant kingdom? Previously, sequence homology-based methods have been widely used for RLK identification and classification, which probably will miss the distantly related proteins but with similar structures and potential novel families unmentioned in the literature. The academic community urgently requires a dedicated database for a systematic overview of the RLK gene family, providing data support for in-depth research on RLK genes. Here, we used a topology-based method to accurately isolate the RLKomes from proteomes. The obtained RLKomes were further classified into (sub)families based on ECD domains. We constructed a comprehensively curated plant RLK database (https://biotec.njau.edu.cn/rlkdb/), which contains valuable resources for investigating the origin and evolution of the RLK family and multiple online tools for personalized analysis.

Auxin and auxin-mediated signaling pathways involved in the regulation of lateral root development are well documented. Although exocytic vesicle trafficking plays an important role in auxin efflux carriers PIN recycling, and polar auxin transport during lateral root formation, however, the mechanistic details of these processes are not well understood. Here, we demonstrate that BYPASS1-LIKE (B1L) regulate lateral root development via exocytic vesicular trafficking-mediated polar auxin transport in Arabidopsis. In b1l mutants, the number of lateral roots increased significantly, and the phenotypes were mainly attributed to lateral root primordium initiation but not to the defects in lateral root primordium development. Furthermore, the auxin signal was stronger in the lateral root primordium of the b1l mutant at stage I than those observed in the wild-type (WT). Moreover, exogenous auxin and auxin transport inhibitory treatments indicated that the phenotype of lateral roots in b1l mutants can be attributed to higher auxin levels and that B1L regulates auxin efflux. Consistently, auxin efflux carriers PIN1-GFP and PIN3-GFP were expressed at higher levels in the lateral root primordium of the b1l mutants. Interestingly, we found that B1L interacted with the exocyst and b1l mutant showed a defect in PIN2 exocytosis. Finally, we found that B1L cooperated with EXO70B1 to regulate lateral root formation. Our findings reveal an essential regulatory mechanism of B1L that interacts with the exocyst to regulate PIN-mediated polar auxin transport and lateral root initiation.

Pyridoxal-5-phosphate (PLP), the enzymatic cofactor form of Vitamin B6 (vitB6), is a versatile compound that has essential roles in metabolism. Cellular PLP homeostasic regulation is currently not well understood. Here we report that in Arabidopsis, biosynthesized PLP is sequestered by specific aminotransferases (ATs), and that the proteins ROOT UV-B SENSITIVE 1 (RUS1) and RUS2 function with ATs to regulate PLP homeostasis. The stunted growth phenotypes of rus1 and rus2 mutants were previously shown to be rescuable by exogenously supplied vitB6. Specific residue changes near the PLP-binding pocket in ASPARTATE AMINOTRANSFERASE2 (ASP2) also rescued rus1 and rus2 phenotypes. In this study, saturated suppressor screens identified 14 additional suppressor of rus (sor) alleles in four aminotransferase genes (ASP1, ASP2, ASP3, or ALANIN AMINOTRANSFERASE1 (AAT1)), which suppressed the rus phenotypes to varying degrees. Each of the sor mutations altered an amino acid in the PLP-binding pocket of the protein, and sor proteins were found to have reduced levels of PLP conjugation. Genetic data revealed that the availability of PLP normally requires both RUS1 and RUS2, and that increasing the number of sor mutants additively enhanced the suppression of rus phenotypes. Biochemical results showed that RUS1 and RUS2 physically interacted with ATs. Our studies suggest a mechanism in which RUS1, RUS2 and specific ATs work together to regulate PLP homeostasis in Arabidopsis.

Brassinosteroids (BR) and Target of Rapamycin Complex (TORC) are two major processes coordinating plant growth and stress responses. BRs function through a signaling pathway to extensively regulate gene expression and TORC is known to regulate translation and autophagy. Recent studies revealed that these two pathways crosstalk, but a system-wide view of their interplay is still missing. Thus, we performed transcriptome, proteome, and phosphoproteome profiling of Arabidopsis mutants with altered levels of either BIN2 or RAPTOR1B, two key players in BR and TORC signaling, respectively. We found that perturbation of BIN2 or RAPTOR1B levels affects a common set of gene-products involved in growth and stress responses. Additionally, we performed Multiplexed Assay for Kinase Specificity (MAKS), which provided a system-wide view of direct BIN2 substrates. Furthermore, phosphoproteomic data was used to reconstruct a kinase-signaling network and to identify novel proteins dependent on BR and/or TORC signaling pathways. Loss of function mutants of many of these proteins led to an altered BR response and/or modulated autophagy activity. Altogether, these results provide genome-wide evidence for crosstalk between BR and TORC signaling and established a kinase signaling network that defines the molecular mechanisms of BR and TORC interactions in the regulation of plant growth/stress balance.

The levels of the important plant growth regulator indole-3-acetic acid (IAA) are tightly controlled within plant tissues to spatiotemporally orchestrate concentration gradients that drive plant growth and development. Metabolic inactivation of bioactive IAA is known to participate in the modulation of IAA maxima and minima. IAA can be irreversibly inactivated by oxidation and conjugation to Aspartate and Glutamate. Usually overlooked because its reversible nature, the most abundant inactive IAA form is the IAA-glucose (IAA-glc) conjugate. Glycosylation of IAA is reported to be carried out by the UDP-glycosyltransferase 84B1 (UGT84B1), while UGT74D1 has been implicated in the glycosylation of the irreversibly formed IAA catabolite oxIAA. Here we demonstrate that both UGT84B1 and UGT74D1 modulate IAA levels throughout plant development by dual IAA and oxIAA glycosylation. Moreover, we identify a novel UGT subfamily whose members modulate IAA homeostasis during skotomorphogenesis by redundantly mediating the glycosylation of oxIAA.

Stomata play critical roles in gas exchange and immunity to pathogens. While many genes regulating early stomatal development up to the production of young guard cells (GCs) have been described in Arabidopsis, much less is known about how young GCs develop into mature functional stomata. Here we performed a maturomics study on stomata, with "maturomics" defined as omics analysis of the maturation process of a tissue or organ. We developed an integrative scheme to analyze three public stomata-related single-cell RNA-seq datasets and identified a list of 586 genes that were specifically up-regulated in all three datasets during stomata maturation and function formation. The list, termed sc_586, is enriched with known regulators of stomatal maturation and functions. We selected two candidate G2-like TFs genes, MYS1 and MYS2, from the list to investigate their roles in stomata. Our results showed that these two genes redundantly regulate the size and hoop rigidity of mature GCs, and their double mutations caused mature GCs to have severe defects in regulating their stomatal apertures. Our analysis thus provides a valuable gene list for studying GC maturation and function formation.

Carotenoids are health-promoting plastidial isoprenoids with essential functions in plants as photoprotectants and photosynthetic pigments in chloroplasts. They also accumulate in specialized plastids named chromoplasts, providing color to non-photosynthetic tissues such as flower petals and ripe fruit. Carotenoid accumulation in chromoplast requires specialized structures and proteins such as fibrillins. Although fibrillins were first reported as structural components of carotenoid sequestering structures in chromoplasts, later work revealed roles in chloroplasts and other plastid types. However, the association of fibrillins with carotenoids in plastids other than chromoplasts has remained unexplored. Here we show that a member of the fibrillin family, FBN6, interacts with phytoene synthase (PSY, the first committed and rate-determining step of the carotenoid pathway) to promote its enzymatic activity. Transient overexpression of FBN6 in Nicotiana benthamiana leaves results in a higher production of phytoene, the product of PSY activity, whereas loss of FBN6 activity in Arabidopsis thaliana mutants dramatically reduces the production of carotenoids during seedling deetiolation and after exposure to high light. Our work hence demonstrates that fibrillins not only promote the accumulation of carotenoids but also their biosynthesis.

Hybridization and polyploidization are pivotal to plant evolution. Genetic crosses between distantly related species rarely occur in nature mainly due to reproductive barriers but how such hurdles can be overcome is largely unknown. xBrassicoraphanus is a fertile intergeneric allopolyploid synthesized between Brassica rapa and Raphanus sativus in the Brassicaceae family. Genomes of B. rapa and R. sativus are diverged enough to suppress synapsis formation between non-homologous progenitor chromosomes during meiosis, and we found that both genomes reside in the single nucleus of xBrassicoraphanus without genome loss or rearrangement. Expressions of syntenic orthologs identified in B. rapa and R. sativus were adjusted to a hybrid nuclear environment of xBrassicoraphanus, which necessitates reconfiguration of transcription network by rewiring cis-trans interactions. B. rapa coding sequences have a higher level of gene-body methylation than R. sativus, and such methylation asymmetry is maintained in xBrassicoraphanus. B. rapa-originated transposable elements were transcriptionally silenced in xBrassicoraphanus, rendered by gain of CHG methylation in trans via small RNAs derived from the same sequences of R. sativus subgenome. Our work proposes that not only transcription compatibility but also a certain extent of genome divergence supports hybrid genome stabilization, which may explain great diversification and expansion of angiosperms during evolution.

Dendrobium officinale is a typical epiphytic orchid. We report the telomere-to-telomere (T2T) genome assembly for D. officinale, representing the first T2T reference genome within the Orchidaceae family. The assembly is anchored to 19 chromosomes and contains 38 complete telomeres and 15 characterized centromeres. We further generated haplotype-resolved assemblies of the autotetraploid genome, identifying 12,761 sets of tetra-allelic genes. Based on synonymous substitution analysis, we inferred that the autotetraploidization event occurred approximately 0.86 million years ago. A systematic analysis of the SWEET gene family across the genus Dendrobium revealed that the gene family size is shaped primarily by epiphytic types and environmental factors. In D. officinale from Langshan, eight SWEET genes were specifically expressed in roots, suggesting they may play specialized roles in the root mycorrhizal system, potentially contributing to the D. officinales ability to recruit and maintain fungal partners. Together, these resources provide valuable foundations for studies of orchid evolution, functional genomics, and molecular breeding.

For the first time, we assembled two highly continuous, completely gap-free reference genomes of the Rubus subgenus: Rubus hirsutus Thunb. Penglei (XMM) and Rubus eustephanos Focke ex Diels Dahongpao (DHP), which are widely distributed in southern China with similar phenotypic traits (Figures 1A, 1B, and Figure S1), yet ripe fruits display distinct sugar accumulation levels (Table S1), making them ideal candidates for investigating the mechanisms underlying sugar accumulation in the Rubus genus. The XMM (213.53 Mb, 28,204 genes) and DHP (218.26 Mb, 28,569 genes) genomes exhibit close evolutionary relationships, diverging approximately 3.21 Mya. Comparative genomics identified extensive synteny, interspecific structural variations (translocations, inversions, segmental duplications), and presence/absence variation (PAV). Using Hi-C interaction heatmaps and Sanger sequencing, we validated interspecific structural inversions. Additionally, we identified a sugar transporter gene (MFS1), which is present in XMM but absent in DHP. Combined analysis of the gene family expansion/contraction and transcriptome identified two conserved key genes (RhSTP13 and RhSTP7) associated with sugar accumulation in Rubus genus and displayed distinct roles through transient expression assay. To facilitate functional genomics study, we also established a comprehensive Rubus database, RubusDB, a freely accessible repository consolidating all genomic, transcriptomic and phenotypic data of Rubus genus. These findings provide a foundational framework for elucidating the genetic basis of sugar accumulation, genome diversification, and trait improvement in Rubus species. O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=170 SRC="FIGDIR/small/646780v2_fig1.gif" ALT="Figure 1"> View larger version (72K): org.highwire.dtl.DTLVardef@1baca20org.highwire.dtl.DTLVardef@2c01eaorg.highwire.dtl.DTLVardef@131dc00org.highwire.dtl.DTLVardef@62a04f_HPS_FORMAT_FIGEXP M_FIG O_FLOATNOFIGURE 1.C_FLOATNO Phenotypic characterization, de novo genome assemblies, variation, and evolution of Rubus hirsutus Thunb. Penglei (XMM) and Rubus eustephanos Focke ex Diels Dahongpao (DHP). (A) Illustration of XMM whole plants, stem, flower bud and flower, respectively. (B) Illustration of DHP whole plants, stem, flower bud and flower, respectively. (C) The chromosome karyotype analysis of XMM. Bar = 10 {micro}m. (D) Circular diagram of XMM and DHP reference genomes. (a) Chromosomes are represented by centromeres (dark red) and telomeres (blue); (b) CEN18 density; (c) CEN17 density; (d) Gene density; (e) TE density; (f) PAV density; (g) SNP density; (h) Collinear lines at the center of diagram highlight homoeologous chromosomes relationships and non-homoeologous regions. (E) Genomic alignments between XMM and DHP. Inversions, duplications and translocations are marked with orange, blue and green ribbons, respectively. (F) Identification of large inversions in chromosome 7 between XMM and DHP. The three heatmaps show the chromatin interaction matrix, including mapping Hi-C data of DHP against DHP genome (left), mapping Hi-C data of DHP against XMM genome (middle) and mapping Hi-C data of XMM against XMM genome (right). The lower panel illustrates gene alignments between XMM and DHP. (G) Synteny, structural rearrangements, monomer, and identity distribution in the 16-20 Mbp region of chromosome 6 between XMM and DHP. (H) Local genome synteny of chromosome 5, with a structural variation that is present only in XMM. The PAV region includes a gene related to sugar transporter (MFS1, XMM05G025540). (I) Gene-level matchings of MFS genes in the PAVs region between XMM and DHP. (J) Phylogenetic relationships and divergence times between raspberries and other Rosaceae species. The black numbers close to the divergence nodes indicate the divergence times and the red bars represent the 95% confidence intervals. The red and black numbers indicate expanded and contracted ortholog groups at the corresponding node. Scale bar corresponds to 10 Mya. C_FIG